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■ Abstract 



Symplectic geometry of the vortex filament in a curved three-manifold is investigated. 
There appears an infinite sequence of constants of motion in involution in the case of 
constant curvature. The Duistermaat-Heckman formula is examined perturbatively for 
the classical partition function in our model and verified up to the 3-loop order. 
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1 Introduction 



Kinematics of a very thin vortex tube in three-dimensional fluid may be described by the 



where 7 = 7 (t, s) denotes the position of the vortex filament in R 3 with t and s being the time 
and the arc-length parameter respectively. 

Hasimoto || introduced a map h : 7 1— > = k exp[i J s r(u)du], in order to transform the 
filament equation into the nonlinear Schrodinger (NLS) equation for ip. Here k and r respec- 
tively denote the curvature and the torsion along 7. Since the integrability of the NLS equation 
was well known, the filament equation was naturally expected to be integrable. Mardson and 
Weinstein Q first described the filament equation as a Hamiltonian equation with the Hamil- 
tonian simply being the length I [7] of the vortex filament. Later Langer and Perline |J used 
this Hamiltonian structure to prove the existence of an infinite sequence of constants of motion 
in involution, and studied the evolution of the vortex filaments in connection with the solitons 
in the NLS equation. 

With this concern in mind, we have investigated the filament equation in a curved three- 
manifold M. Although Langer and Perline have limited M to R 3 , we find an analogous in- 
tegrable hierarchy in the case of constant curvature. We further study the classical partition 
function for the vortex filaments 



space T is neither finite dimensional nor compact, and furthermore because the Hamiltonian 
flow may not be periodic. But the perturbative calculation in our mode reveals that the loop 
corrections to the formula vanish up to the 3-loop. 

2 Integrability 

We begin this section by describing a symplectic structure for the vortex filament in a three- 
manifold M equipped with a Riemann metric g. Everything is considered in the smooth 
category for simplicity. Let T be the space of vortex filaments with fixed end points p and 
q; T is the quotient space of {7 : [0, 1] — > M | 7(0) = p, 7(1) = q} with the reparametrization 
of 7. Hereafter 7 denotes the representative for which the parameter x G [0, 1] is a multiple of 
the arc-length s, namely 



filament equation in the local induction approximation [1], ■ It is formulated as 

c>7 d'y <9 2 7 



(1) 



dt ds ds 2 ' 





(3) 
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is independent of x. Here ( , ) denotes the inner product on the tangent space T^mM. One 
can identify the tangent space T 7 T with the subspace of T(7*TM), and expand X G T(j*TM) 
in the Frenet-Serret frame along 7 such that 

X = fT + gN + hB, (4) 

where T is the unit tangent vector to 7, N is the unit normal vector and B is the unit binormal 
vector. Let I [7] be the length of 7, so that s — £ [7] x. The Frenet-Serret equations are 

V s T = kN, V s N = -kT + tB, V s B = -tN, (5) 

with V being the connection on 7*TM induced by the Levi-Civita connection on TM. Let p 
be the projection from T(7*TM) to T 7 T, then one can show that the tangent component of 
v = p (X) G T 7 r satisfies 

* VT = i - 1 (V m v,p-)+£KV H , (6) 

and (V x v , d^y / dx) is a constant. Fixing this constant by the boundary conditions X(0) = 
X(l) = 0, one obtains 

p(X)—v — £( nv^dx — x / KVfidx) T + t^N + t>eB . (7) 
Jo Jo 

Geometrical structures on T were first studied by Marsden and Weinstein 0] for the vortex 
filament in R 3 , and generalized to the loop space for a three-manifold M by Brylinski [^]. It is 
straightforward to find those for the vortex filament in M. 

i) Complex structure 

For the tangent vector v G T 7 r, J generates the 90-degree rotation 

J(v) = -p(Txv) , J 2 = -l. (8) 

Choosing (t> N , t> B ) as coordinates for T 7 r, we find that J corresponds to the multiplication by 
i for the complex function ^n(^) + ivb(x). Hence J induces a complex structure on T. 

ii) Riemann structure 

The Riemann structure on T is simply defined by 

(u,v)r = £ f (u n vn + u b vb) dx (9) 

J 

for T^T, and satisfies the hermitian condition 

(u,v) r = (J(u),J(v)) r . (10) 
Note that even though ( , )r ignores the T-components, it is non-degenerate. 
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iii) Symplectic structure 

The volume form v on M associated with the Riemann metric g provides the symplectic 
structure on T, namely 



u(u, v) 



o 



v(-r~, u, v) dx. 



dx ' 



Using the Frenet-Serret frame, one can rewrite this as 

u(u,v)—£ [ (u N v B - u B v N ) dx, 
Jo 

which is equivalent to the one constructed from the above two structures 

u(u,v) = (u,J(v)) r . 



(11) 



(12) 



(13) 



Having set out the basic structures, we now turn to the Hamiltonian flows for the vortex 
filament. Let I : T > R be a smooth Hamiltonian function, then the Hamiltonian vector field 
Xg has the form 

X, = J(grad£). (14) 
Choosing ix e u = dl and putting v = dj t /dt \ t =o, we get 



v£[y] 



d 



_ (<ht <ht\ 
dt Jo \ \ dx ' dx J 

[l] fo ' dx) 



dx 



t=o 



dx, 



(15) 



£[y] / (v,Kti)dx. 



grad£ = — p(n N) follows, and therefore 



X f = k B. 



This yields a natural generalization of the filament equation in M 

dt ox ox 

The evolution equations for k and r are the followings 

3k 

dr 
di 



k Ric (B, N) 



(2r— + k— 

dx dx' 



r Ric (T, N) + r 1 J^/, 2 + r 2 K~^ - r 2 + p(T, B)), 



(16) 

(17) 

(18) 
(19) 



where Ric and p denote the Ricci tensor and the sectional curvature on M respectively. In the 
case of constant curvature, these equations take simpler forms 

dK 

dt 
dr 

dt 



^ dx ^ dx ' 



(20) 
(21) 



which turn out to be identical with the equations appeared in || . Hence we can get the following 
proposition for a three-manifold with constant curvature. 

Proposition 

(a) The filament equation is transformed into the NLS equation by the Hasimoto map. 

(b) There is an infinite sequence of constants of motion. 

(c) These constants are in involution. 

Proof 

We assume that k, t and their derivatives of arbitrary order vanish at the boundaries. 



Then it is straightforward to prove (a) and (b) due to the evolution equations (p0| ) and (21). 
Using the explicit form of the Hamiltonian vector fields X n (see Remark (2)), we can confirm 
the commutativity ui(X n , X m ) = for any n and m with the help of |J, and consequently prove 
(c). 

Remarks 

(1) The constants of motion are as follows |J: 

J_ 2 [ 7 ] = £fr], I-i[y]=e f'rdx, (22) 

Jo 

J n [ 7 ] = J„o%] (n = 0,1,2,...), 

where h is the Hasimoto map h[y] = k exp[i£ [7] Jq rdx], and J n 's are the constants of 
motion in the NLS equation pi given by 



lM=t t\^J n {^,i>)dx, (23) 
Jo 2 

and 

d ~ 1 - n ~ 

J = ij, J n+1 = -i —J n - -tpY^Jk-iJn-k- (24) 

(2) We find mutually commuting Hamiltonian vector fields X n for I n [y]: 

X_ 2 = kB, X_i = i2X_ 2 , (25) 
X n = R n+2 X_ 2 -cR n X^ 2 (n = 0,l,2,...), 

where c denotes the constant curvature and R the "recursion operator" defined by 

R(v) = -f x p(T x V x v) (26) 

for v G T 7 T. R coincides with the one appeared in || when we restrict v to the Hamilto- 
nian vector fields. 
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(3) Langer and Perline interpreted the Hasimoto map as a Poisson map between the Poisson 
structure on the space of the vertex filaments and the "forth" Poisson structure on the 
space of the NLS fields. We have found no such correspondence in our model, because the 
deformation of the vortex filament also changes its length £. In the case of || , however, 
the vortex filament extends boundlessly, so that the arc-length parameter is simply a 
parameter and does not change under the deformation. A different approach to the 
integrability of the vortex filaments has been investigated in [0] recently. 

The filament equation belongs to an infinite hierarchy of Hamiltonian systems {d'y/dtn = X n | 
n = —2,— 1,0,...}, and all Hamiltonian flows in this hierarchy are transformed into those in 
the NLS hierarchy. In fact, the differential of h yields 

dh : X n i — > X n+A - 2cX n+2 + c 2 X n (mod iip), (27) 

where X_2 = X-i = 0, and X n (n = 0,1,2,.. .) are the Hamiltonian vector fields associated 
with I n [i>), i.e., X n = — zgrad/„; first two are 

dh{X. 3 ) = z(^ + i|^| 2 ^), (28) 

3 Classical partition function 

In this section we evaluate the classical partition function (§) with T> / y being the symplectic 
volume form on T. The stationary phase method provides an asymptotic expansion for Z{0) 
as (3 I— > oo, such that 

Z((3)= E Z WKB [ 1 ,^(l + ^- + ^ + ...). (30) 

grad^[ 7 ]=0 P P 

The exactness of the stationary phase (WKB) approximation has been of interest due to the 
Duistermaat-Heckman formula pi, where they have shown that if T is a compact symplectic 
manifold and I is a periodic Hamiltonian with isolated critical points, WKB approximation 
becomes exact for (El), i.e., the asymptotic expansion terminates at Z WKB . In more general 



arguments presented in ||TT[, the fixed points are not necessarily isolated, and it is not manda- 
tory to consider the circle action alone according to the analogous results obtained for higher 
dimensional tori. For the infinite dimensional symplectic manifolds, the WKB exactness has 
not been proved rigorously, but a "proper" version of WKB approximation should yield a re- 
liable result for a large class of integrable models [12, [l^, 14, Tf|. With this notion in mind, 



we present the explicit calculation of the asymptotic expansion (p0|) . For simplicity, we will 
assume the followings: 
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(1) M is a three- manifold with a constant curvature c, so that the filament equation is 
integrable in the sense of Proposition. 



(2) Two points p and q on M are not conjugate. Consequently, the Hamiltonian I is a Morse 
function on T, i.e., critical points are the geodesies on M connecting p and q, and further 
the Hessian operator H 1 at each geodesic 7 is a non-degenerate Jacobi operator 



# 7 = -v.v, - dYif 



(31) 



Let us first expand the Hamiltonian I around a geodesic 7. As we can see in ([5]), the cur- 
vature along the geodesic vanishes identically, and £(x) G T y ^M thus satisfies the condi- 
tion (£, T) = 0. Using an infinitesimal deformation of 7 generated by the exponential map 
7 S (x) = exp^r x \[s £ £(x)], we can find the expansion 



00 s n d n 



WE 



00 ^2n r \ 



„=0 W J ° 

Here the integrand W 2n is given by the Bell Polynomial Y r 



s=0 

W 2n (0dx. 



with 



W 2n (0=Y 2n (f 1 ,...,f 2n ,g 1 (0, 



. 1 (2m-3)!! 



pj] , namely 



/m ( 

<7 2 (£) = 2(V^,V x O-2c£ 2 (e,0, 

^(o = (-) m 2 2m - i {(^ 2 ) m (e,o m 

+ (c^^ce.o^Kv.e.o'-ce.ocv^v.oi} (™>2), 

SWl(6) = °- 

First few are given by 

W = l, W 2 = f l92 , 

W 4 = fx g 4 + 3 f 2 gl W 6 = hg 6 + 15 f 2 g 4 g 2 + 15 f 3 g%. 

Now let us evaluate the WKB partition function 



tw 2 {i)dx = (£,# 7 (0>r- 
Jo 
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P£[i\ 



W 2 {0dx 



(32) 
(33) 

(34) 



(35) 



(36) 



(37) 
(38) 



Using the zeta-function regularization technique, we can perform the infinite dimensional inte- 



gral in fl37D , and obtain fl7 



[ 7; /3] = e -/3%]±i*(to(o)-Cff(o)) e iC^(o) (/3^[7]) - ^(°), 



(39) 



with t]h(z) and Cff(^) I- 2 e Q being eta and zeta functions associated with the Hessian operator 
if 7 respectively. Evaluating these functions for 7 with the Morse index ^(7), we find 



W (0) = -l-2//(7), Ch(0) = -1, Ch(0) = 21u 



2sin(v^£[7]) 



and eventually this gives us an explicit expression 
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Z WKB [7,f3] = -e 



1 P -Ptb] 



g=Ff«M(7)_ 



(40) 



(41) 



Since ^(7) is an even integer, the last factor contains no ambiguities. 

We now proceed to the higher-order calculation. It is convenient to choose an orthogonal 
frame {ei, e 2 } along 7 such that 



V x e i = 0, (T, ei ) = fori = 1,2. 



(42) 



In this frame, the kernel of the Jacobi operator if 7 becomes diagonal, and both of the diagonal 
elements are identical to the Dirichlet Green function 



G{x,x') = 2]T 



70 sinfnTrx) sinfwrrr') 



n=l 



(nir) 2 — X 



(43) 



with A = c£ 2 . The 2-loop amplitude a\ = —(3 2 (W4/4!) consists of four diagrams depicted in 
Fig. 1, and those are respectively 



A / G(x 
Jo 



,2 = ^_3 v / Xx + 3 Ax2 



(b) 



(d) 



X £ G(x)G"(x) = -±- 1 -VxX+ 1 -XX' 

1 G"(x) 2 = X 2 f l G(x) 2 , 
J 



where X = coty/X, G(x) = G(x,x') \ x=x >, G'(x) = (d/dx) G(x, x') \ x=x > and G"(x) = 
(d 2 / dx dx') G(x, x') \ x=x i. While G(x) and G'{x) are convergent, G"{x) diverges at the bound- 
aries, thus we have found (c) and (d) by executing the x-integration first and then by regularizing 



S 



the infinite n-summation in terms of the following analytic continuations: 



l .-2» , F ( S ~ D -2s+l 



£(n 2 + a 2 )* 2° 2 I/O) 



a 



i 



4 ^(2.™,), (44) 

_ ^ r ( S - |) -2s+3 



E ( n 2 + a 2), 4 r(s ) 



a 



7T2 9 /7rn\ s 2 



r(s) v a y 

~ n 4 3 ir(s-f) _ s+ b 
> = — 7r 2 — a 2 

^( n 2 + a 2 )s g r(s ) 



E- 2 - *.-*(2™a), 



3 7t5 ~ /7rn\ ; 



4 ifSrXHTj ' ^t<^) 



n=l 

7T 2 



n=l 



- s S?»'(7l K -^ n "> (46) 



4 ^f-«(=) H MP«-). 

where K u (z) is the modified Bessel function. Multiplying (a) through (d) with the weights of 
the diagrams, we conclude that the 2-loop amplitude vanishes. Beyond the 2-loop, however, 
we ought to generalize the analytic continuation for a multiple infinite summation. One might 
think that applying the analytic continuation method directly to the Green function, we could 
regularize the Green function, and thereby making all loop amplitude finite. This is certainly 
true, but regularizing the Green function in this way, we also eliminate the necessarily singu- 
larity at x — x', and obtain non- vanishing 2-loop amplitude as a result. We may avoid this 
difficulty by treating G(x,x') as a distribution w.r.t. x. Let us first examine this on the 2-loop 
and check if the amplitude vanishes. Since G(x, x') may naturally be extended periodically 
(period 2) to R as a function of x, one can redefine it as a distribution G(x, x') such that 

G(x, x') = p-i — = V {sia.[y/\ (x - 2n)\ (x' - 1)] H(x; 2n, x' + 2n) 

V A sin V A neZ L 

+ sinf-v/Aa:'] sinfv^ (x - 2n - 1)] H(x; x + 2n, 2n + 1)} (47) 



+ ty 1 7t E { x } > 

VAsm VA neZ 
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where H(x; a, b) denotes the characteristic function for the interval [a, b] C R. Similarly G(x) 
may also be extended periodically (period 1) to R 

G(x) = j=^- — = {sintv 7 ! (x - n)] sinfv 7 ! (x - n - 1)] H(x; n, n + 1)} . (48) 

VAsin vA neZ L J 

Using the periodic delta function 5(x; n) (n is the period), we may evaluate the second derivative 

— — G(x, x') = ^—= Y f cosh/A (x - 2n)l cosfVX fx' - 1)] #(x; 2n, x' + 2n) 

dzdx' sinv^,^ 1 J 

+ cos[VAx'] cos^ (x - 2n - 1)] if(z; x' + 2n, 2n + 1)} (49) 

VA 



and similarly 



V { x -> -x } + 5(x - x'; 2) + 5(x + x'; 2) 
sin V A n6Z 



G"(x) = -A G{x) - V\ cot yf\ + i 5(x; 1). (50) 



The delta function appears only in G"(x), and we confirm the vanishing of the 2-loop amplitude 
by using 

/ dx5(x;l) 2 = 5(0;1) = 0, (51) 
Jo 

which is consistent with the ^-function regularization because of 5(0; 1) = 1 + 2((—l). 

The 3-loop amplitude a 2 = j3 3 (j3(W 4: /4:\) 2 /2 — W e /6\) consists of 30 diagrams depicted 
in Fig. 2. Evaluating them by means of (fj?D, (13) and 5(0; 1) = 0, we find that 29 diagrams 
contain no ambiguities due to the the integration formulae 

fdx [ X dy5(x; l)5(y; 1) F(x,y) = i F(0, 0) + - F(l, 0) + ~ F(l, 1), (52) 
jo jo 8 4 8 

/' rfx f dy [5(x - y; 2) + 5(x + y; 2)] [5(x; 1) + <%; 1)] F(x, y) 
jo jo 

= If(0,0) + ^F(1,1), (53) 

f'dx f X dy[5(x-y;2) + 5(x + y;2)] 2 F(x,y) = \ F(0, 0) + J l). (54) 
Jo jo 8 8 

Here the last equality follows from 5(0; 2) = 0. Yet, in the diagram whose weight is —480, we 
encounter an ambiguous integral 

f dx f dy [5(x -y;2) + 5(x + y; 2)]5(x; 1) F(x, y) = p F(0, 0) + q F(l, 1), (55) 
JO jo 

where p + g=l/2asis shown in (53), but p or q alone cannot be determined unless we specify 
the regularization of the delta function. If we were able to define the analytic continuation of 
the infinite double sum, this ambiguity would not appear, but we have no choice at our hand 
other than putting q = 1/16, and obtain the vanishing 3-loop amplitude as a result. 
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Ambiguities appearing in higher loops are inevitable, because they relate to the regular- 
ization ambiguity of the integration measure X>7, which has never been defined rigorously in 
the first place. Both methods we have presented here reveal that the degree of ambiguity gets 
larger as the order of loops increases. In the analytic continuation method, ambiguity arises 
from the variety of the analytic continuation applicable to the multiple infinite summation, 
whereas in the distribution method, the delta-function integration, particularly the finite part 
of the boundary contribution, is the source of the ambiguity. Nevertheless our lower order 
calculations suggest that by regularizing order by order, one can eliminate all higher loop 
corrections, and thereby preserving the Duistermaat-Heckman formula. The symplectic struc- 
ture has been studied thoroughly in compact finite dimensional manifolds, but little is known 
for the infinite dimensional ones, which include most of the integrable hierarchies. This is ex- 
actly the place where the physical interests are, and the Duistermaat-Heckman formula would 
throw a new light over the integrable hierarchies as we have caught a glimpse of it here. 

The authors would like to thank Dr. N. Sasaki for helpful discussions. 
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Figure captions 



Figure 1. 2-loop diagrams. "Dot" denotes a derivative on the Green function; for instance, 
the propagator with one dot represents G'(x), and the one with two dots G"(x). The attached 
numbers are the weights of the diagrams. 

Figure 2. 3-loop diagrams. Since there appear no second derivatives in (|35| ) and the number of 
derivatives is always even, double dots on a single propagator must go to the separate vertices 
and the number of the derivatives at each vertex must be even. One must interpret dots 
accordingly for the diagrams with two vertices. Note that for a couple of diagrams in the 
top group and for a couple in the middle, though the resulting diagrams are inequivalent, this 
simple rule does not tell to which vertex dots are supposed to go. In those diagrams, dots are 
placed closer to the vertices to which they are supposed to go. 
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